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Abstract 

The U (N) lattice gauge theory in 2-dimensions can be considered as the 
statistical mechanics of a Coulomb gas on a circle in a constant electric 
field. The large limit of this system is discussed and compared with ex- 
act answers for finite N. Near the fixed points of the renormalization group 
and especially in the critical region where one can define a continuum the- 
ory, computations in the thermodynamic limit (N — )■ oo) are in remarkable 
agreement with those for finite and small N. However, in the intermediate 
coupling region the thermodynamic computation, unlike the one for finite A^, 
shows a continuous phase transition. This transition seems to be a pathol- 
ogy of the infinite limit and in this simple model has no bearing on the 
physical continuum limit. 



*Work supported in part by the NSF: Grant No. PHY- 78-01224. 
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Introduction 

The 1/A^ expansion has proved very useful in studying the physical spec- 
trum of N component spin systems (cr-models) and several other tractable 
field theories with large dimension of internal symmetry group, in two dim.^ 
It is used on the assumption that the qualitative, and, in certain cases, even 
quantitative character of the spectrum is independent. Then taking the 
large N limit enormously simplifies the computation since one can use a 
combination of a steepest descent and mean field approach in which each 
component of the spin interacts with the mean field of the other compo- 
nents. 

For SU{N) color gauge theories of the strong interaction, the ex- 
pansion was introduced by 't Hooft^ who showed that in the Feynman graph 
expansion of the theory to each order in g'^N (fixed as N ^ oo), the lead- 
ing contribution comes from planar diagrams. Diagrams with holes (quark 
loops) and handles are suppressed by factors of and respectively. 
In this sense, the expansion in is analogous to the topological expansion 
of the dual model and as ^ oo the mesonic and gluonic states appear 
as free particle states. Witten^ has recently incorporated baryons as soli- 
tons into this scheme and has enumerated and emphasized the qualitative 
agreement between simple and basic aspects of hadronic spectroscopy and 
the 1/N expansion. However, the computational difficulties of planar QCD 
remain intractable in the diagrammatic language of its formulation, though 
some remarkable progress in the counting of planar diagrams for quartic and 
cubic vertices has been made by Brezin et al., using a WKB approach in 

Now in the statistical mechanics of component spin systems the rela- 
tion of the large N limit to a mean field computation has received precise 
mathematical formulation by the introduction of a Lagrange multiplier as 
an auxiliary field. It is desirable to have an analogous formulation in the 
case of QCD, which avoids direct encounters with the diagram technique. 
Further, it is desirable to evaluate the relevance of the large N limit. The 
question is very simple: to what extent and in which domains of coupling 
is a thermodynamic computation in which N ^ oo, relevant for the case of 
finite N. 

We present in this paper a study of U{N) lattice gauge theories in 2- 
dimensions to clarify these questions. Our method follows studies in the 
statistical theories of spectra and ref. (4). 
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In two dimensions, since plaquette variables are the independent ones, 
the partition function reduces to that of a random unitary matrix with N"^ 
degrees of freedom. This in turn is the partition function of a Coulomb gas on 
a circle in a constant electric field E = 2(3 /N . The problem of computing the 
free energy and other thermodynamic quantities as N ^ oo reduces to the 
solution of a singular integral equation for the charge density, with a Cauchy 
kernel. This problem is always exactly soluble. The solution depends on the 
strength of the electric field E. For < < 1 there is no gap in the charge 
distribution and the charge density is conjugate to the electric field. For 
1 < < oo the charge distribution has a gap and the problem reduces to 
the Riemann-Hilbert problem for a simple arc in the plane. The appearance 
of a gap in the charge density signals a continuous phase transition sX E = 
1. An order parameter that measures the randomness of the system has 
a constant value in the 'strong coupling' phase < ii^ < 1 and goes to 
zero for large E in the spin wave phase. The string constant is computed 
and a simple renormalization group used to discuss the physically relevant 
continuum limit. 

The Coulomb gas partition function is also a Toeplitz determinant in- 
volving Bessel functions. For finite N it is obviously an analytic function of 
the temperature ^. Taylor expansions around the fixed points ^ = and 
(3 — oo indicate an excellent agreement with the computation for — > oo. 
To us, this agreement is far from self-evident. Near (3 — oo, this can be seen 
to result from a scaling argument. 

An interesting mathematical corollary is that the problem of evaluating 
Toeplitz determinant of order N, can always be mapped to the problem of 
a Coulomb gas on a circle, which in the N — oo limit is always exactly 
soluble. Toeplitz determinants are encountered in the Ising model and their 
evaluation is usually done by using Szego's theorem (c.f. ref. (11)). 



4 



I. Equivalence of U{N) Gauge Theory in 2-dimensions to the Sta- 
tistical Mechanics of a 1-dimensional Coulomb gas 



The partition function of the U (N) lattice gauge theory is defined by^ 



Z(V,N,I3) = J d/iexp 



f3j:{U{P)-N) 



(1) 



dfi is a real measure over the gauge group elements at each link; the summa- 
tion over P is over all oriented plaquettes and U (P) is a product of group 
elements, in the fundamental representation, at the 4 oriented links of P. 
U'^{P) would be the group element associated with the finks in the opposite 

direction. /3 = — 5- is the temperature and V is the volume of the system in 

lattice units. 

In two space-time dimensions it is possible to treat plaquette variables 
as independent.^ A simple way to see this is to fix the axial gauge in Z: all 
time-like links are chosen to be one and all space-like links at one particular 
time 'to' are chosen to be one. Since the time links are trivial the space- like 
links at time can be expressed as a unique product of a time ordered string 
of independent plaquette variables, the expression for Z in (1) becomes an 
integral over a single unitary matrix 



V 



g -2/3ivy 



Hence, the problem is reduced to the computation of the partition func- 
tion of a random unitary matrix with A^^ real independent variables 

z = j dC/e^(^^+^^+) (2) 

Integration over the group U{N) is known. It is convenient to make a sep- 
aration of variables into group invariant and non-invariant parameters. Let 
R be the unitary matrix that brings U to diagonal form 

U = R+DR, Dij = Sije'^i 

R is arbitrary up to an invariant subgroup of U{N), hence, after appropriate 
choice of gauge it depends on iV^ — N parameters and 

N JO ^ 9 

dU = CNdR n — n - e''^' (3) 



=1 27r 



i<3 
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Substituting (3) in (2) we get 



1 

m 



2- Jt dOi 

IT— e 



iV2 



^1 1 

fT.j;^^os9, + 2j: — log\e^^ 



i<3 



7V2 



(4) 



The factor N\ appeared by choice of Cn in (3) such that z 
Now define an action 



1 for ^ = 0. 



3 ^ 1 



NUN 



Kj 



7V2 



then 



1 

AH JO 



2^ de, 

H 27r 



log 



(5) 



(6) 



The action 5" = o(l) in N provided ^/N — E/2 is a fixed parameter of 
the theory. If so, then for large N we can consider evaluating 2; by a WKB 
expansion in the small parameter 1/iV^. 

The action S in (5) is essentially the electrostatic potential energy of 
N charges of strength 1/N distributed on a unit circle in the plane, in the 
presence of a constant electric field E/2 along the x-axis. The repulsive 
logarithmic interaction tends to distribute the charge uniformly about the 
circle and the electric field tends to drive all — > 0. There is another 
analogy with spin systems. Define a unit spin vector 



(Tj = (cos 6i, sin 9i) 



then 



H 
~2 



N 



i=l 



1 



-S ^ ^.af + y — \og{l 



^ iV2 



^3) 



which is the potential energy of unit spins equally spaced along the x-axis in 
a constant magnetic field H/2 — /3/N. The logarithm represents a con- 
tinuous version of an anti-ferromagnetic coupling. In the absence of H, 
(an) = (cos sin In the rest of this paper we shall only consider 

the electrostatic analogy. 
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II. WKB and Mean Field Calculation 



If we define the free energy per degree of freedom f{E) hj z = exp{N^f{E)), 
then we have the WKB expansion in 

/(B) = -25„ + §+0(ij) (7) 

5*0 is the minimum value of the potential energy S and 5*1 is the result of 
gaussian fluctuations around the most probable configuration that minimizes 
S. For the present, we shall be concerned with evaluating F{E) as N ^ oo. 
As we have noted, the point charges on the unit circle carry a charge 

— 1/N. The interaction energy between two charges is (^jp^ ; which means 

that they are very weakly coupled. However, the interaction energy of any 
one charge in the mean field of the other — N — 1 iv charges is 
0(1). This motivates a mean field theory type computation''' in which one 
introduces a macroscopic charge density u{9) > 0, on the unit circle, such 
that u{9)d9/2n is the fraction of the number of charges between 9 and 9 + d9. 
By definition 

d9 

L («' 

In this continuum limit as N ^ oo we have 

/? /.27r rlf) n2TT rig n2n rlA 

-S[u] = I / ^u{9) COS9 + P ^ %{9) log \e^' - e^^|«(0) (9) 

Z Jo ZTT Jo ZTT Jo ZTT 

P stands for the Cauchy principal value. 

The most probable charge distribution is a stationary point of (9) sub- 
ject to the normalization (8). Hence, this distribution satisfies the following 
singular integral equation 

E /"St rl9 

^cos(l) + P —u{9) \og\e'^ -e"t'\^X 
2 Jo 271 

A is a lagrange multiplier corresponding to (8). It represents the average 
potential at due to all the other charges, and it is 4> independent. Taking 
(p derivatives we have 

£;sin0 = -P l""" ^u{9) cot (10) 
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which states that the tangential forces on the charge at due to all the other 
charges balance the force of the electric field. 

III. Solution of Integral Equation 

Let us begin by assuming that the function u{9) is smooth in the sense 
that u'{9) is continuous for < 6* < 27i. The solution to (10) can be im- 
mediately written down if we recall the Hilbert transform formula for the 
circle: 

If ^{Z) is analytic for \Z\ < 1, then on the unit circle we have 



The choice ^{Z) = 1 + EZ gives the correctly normalized solution 



The condition u{9) > for < ^ < 27r, however, requires < < 1. 

When E = 0, u{9) = 1 as expected; turning on the electric field along the 
x-direction diminishes the density around 9 = tt and produces a crouching 
of charge near ^ = 0. Also u{9) = u{—9). Note that u{9) is smooth and 
never vanishes for £■ < 1. At = 1 the distribution hits a zero at the single 
point 9 = 71. This gives us a hint that as E starts exceeding the critical 
value 1, there would start appearing a gap in the charge distribution around 
6* = TT which would increase as a continuous and monotonic function of E, 
1 < E < oo. In the following we shall estabhsh this picture by solving the 
integral equation (10) with the boundary condition that u{9) — outside the 
interval —a < 9 < a. For a 7^ tt we do not assume u{9) to be smooth at the 
end points 9 — a. The integral equation now becomes 



To solve (12) we introduce complex notation t — e'^ and to — e"^. (12) 
becomes 




u{9) = 1 + E cos 9 



(11) 





(12) 
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^i;iu(t> = l, U=ti = e-" (13) 

Now the important point is that the normahzation condition enables us to 
write the integral equation in the form 

P /■tz dt 

TTl Jti t — to 

The kernel is Cauchy type and we have adapted the general formalism to solve 
such equations to our problem.^ One starts out by considering the function 

analytic in the plane cut by the (ti,t2)- The direction of the arc is as- 
sumed counterclockwise. Denote the boundary values of this function as z 
approaches a point t on the arc {t^ti, t2), from the left and right by ^'^{t) 
and $~(t) respectively. Then by the Plemelj formulae^ for the boundary 
values of analytic functions we have 

- <l>"(t) = 2M(t) (15) 

$+(t) + ^-{t) = - r f = E{t - t-') (16) 

Equation (16) specifies the Reimann-Hilbert problem for the arc (^1,^2): to 
find an analytic function in the cut plane whose boundary values satisfy (16). 
In our case we have an additional condition to satisfy at 2; = 00 which follows 
from the normalization of u{t), 

$(00) = -1 (17) 

The discontinuity formula (15) (a refiection of Gauss' law in electro-statics) 
implies that a solution to the Reimann-Hilbert problem solves the integral 
equation (13). 

The function h{z) = [{z — ti){z — ^2)]^''^ clearly solves the homogeneous 
problem 0+(i) -|-0'~(i) = up to an entire function. This choice of the square 
root function vanishing at the end points leads to u{ti) — u{t2) — 0. Now 
write 

^z) = h{z)H{z) 
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substitution in (16) gives 

H+{t)-H-{t) 



^J{t-h){t-t2) 



(18) 



(By the square root in the discontinuity formula for H[z) we will mean 

1/2 

.) Prom (17) we know that H{oo) ~ 0. 
;he Plemelj formulae 

E{t-t-^) dt 



h+(t) = expi6/2 sin^ - - sin^ - 
Hence H{z) is uniquely given by 



H{z) = / 



This integral can be evaluated by standard methods. Note 



dt 



-^r ^(t-t,){t-t2)-t-^ 

where F encloses the cut (^1,^2) and z lies outside it. This is always possible 
since z / (^1,^2)- H{z) has poles at t = and t = z with residues E/2z and 
E{z — z~^)/2[{z — ti){z — t2)Y^'^ respectively The contribution from infinity 
is E/2 where we have taken the branch 



Hence 



and 



\J{z — ti){z — t-i) —z as 2; — >■ 00 
E{z) 



E 
~2 



I =+1+0 
^^^_t,){z-t2) 2 



$(00) = —1 implies 



1 — cos a 



1 

'e 



a 



sm 



2 E 2' 

Prom the discontinuity formula (15) we get 

-{-■ 

2 \t 
10 



i< £; < 00 



u{t) 



9 
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u{9) = 2E cos -\ E-^ - sin^ - > 



We summarize the computation 
' l + £;cos^ 

u{e) = < 



2£;cos-V£^"^ - sin' 9/2 



0<E<1 
KE <oo 



(19a) 
(1%) 



Our expectation for the appearance of a gap at 6' = tt beyond ii^ = 1 is 
confirmed; also the size of the gap 2 a is a continuous and monotonic function 
ofE: 1/E 

IV. Computation of Thermodynamic Quanti- 
ties in the N ^ oo Limit 

The appearance of the gap at = 1 is reflected in the behavior of the 
thermodynamic functions. From (7) and (9) we have for the free energy per 
degree of freedom 



f{E) 



E^ 
4 

E-^\ogE-^ 



0<E <1 
KE <oo 



(20a) 
(206) 



f{E) and its first two derivatives arc continuous ai E = 1, indicating a 
continuous phase transition at the point E — 1. 

A useful order parameter^ that characterizes the phases is 



N 



R^Y — sin^ 9i 



N 



-Tr 



U-U+ 
2i 



(21) 



In the continuum approximation as A'^ ^ oo we have 

1 ^ /-St H9 

R{E) = lim (- E sin^ 9,) = / -u{9) sin' 9 



Therefore 



f 1/2 



RiE) 



I E 2E'^ 
11 



< £; < 1 

i< £; < oo 



(22) 



(23) 



R represents the average value of the sum of the squares of the y-co-ordinates 
of the charges. R{E) and R'{E) are continuous dX E = 1. In a sense, i? is a 
disorder parameter because it is a measure of the randomness of the system. 
In the phase < < 1 where its expectation value is a constant, the 
distribution of eigenvalues of the unitary matrix is truly random and the 
eigenvalues go over their entire range. We shall call this phase the random 
phase (RP). In the phase 1 < < oo we see that for ^ 1, {R) ^ 0; the 
random distribution has continuously disappeared and the eigenvalues are 
confined to a small region around U = 1. Also the eigenvalue density (196) 
for £■ ^ 1 is written in scaled form 



„W«2v^\/l-<^. (24) 



We call this phase the familiar spin wave phase (SWP). We are aware that 
for £■ > 1, the term SWP is certainly not very appropriate. 

V. The Wilson Loop 

An important order parameter of the original U {N) theory is the Wilson 
loop for a closed curve C, 

W[C] = ^ / ci/zTr (j[ U,^ e^i-^™^-^^ (25) 

Once more, because in two dimensions plaquette variables are the indepen- 
dent variables, we have for a simple curve 

W{C) = 00^ 

uj^^f ci[/(Tr[/)e'^(^^+^^^) (26) 

A is the area enclosed by the curve in lattice units. In the diagonal repre- 
sentation 



1 ^ 



^Ecos^.^ (27) 



which in the continuum approximation of the large N limit becomes 



lim (— y^cos^j)= / — m(6')cos^. 
^-^~\A^^ / Jo 27r ^ ^ 
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Using (19) we have 

r(f)^ o<E<i 

Wa=\ ^ 

i (l - 2^) l<E<OC 

The string constant in lattice units is defined by 



(28) 



1 



WA = e ^""0 



1 



47rlogE/2 
1 1 

^ V 2E 



0<E<1 
1< E <oo 



(29) 



it is the exact analogue of the correlation 'length' in spin systems. Up until 
now we have worked from the statistical mechanics point of view in which 
the basic unit of length is the lattice spacing 'a'. In order to make contact 
with field theory^*^ where the basic unit of length is usually set by the inverse 
renormalized mass we go over to laboratory units in which the string constant 
has the dimension of area 



a{a; E) = a'^ao{E) 



(30) 



This physical quantity is a renormalization group invariant i.e., if a ^ Xa, 
A > 

a{Xa; E) = a{a; Rx{E)) (31) 
Rx{E) is the renormalized coupling. Then from (29) we have 



Rx{E) 



1 



1 Xl/A2 

^~ 2E. 



0<E <1 
l<E<oo 

1 



(32) 



Now consider moving towards short distances, e.g., A = - < 1. The recursion 
formulae which follows from (32) have two trivial fixed points E* — and 
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E* = oo and the flow is away from E* ~ towards E* — oo. Note that 
E — 1 is not a fixed point. At the fixed points we have 

ao{E* = 0) = 
ao{E* = oo) = oo 

Hence £■* = oo is the critical point of the statistical mechanics problem 

for fixed a. In the neighborhood of the critical point E* = oo, the lattice 
constant is an irrelevant parameter and it is possible to define a euclidean 
invariant theory characterized by the adjustable parameter 

aR = a{a -^0,E ^ E*) 

From (29) we have for near the critical point 

Ea^ _ Er 



aR 



27r 27r 



In the field theory notation — = ^ = ~277' hence 



E _ 13 _ 1 



1 

OiR 



nN{gi/a') 

The renormalized parameter aR ranges from to oo, in contrast to the bare 
parameter ckq which ranges from l/27r to oo. This agrees with 't Hooft's 
computation^ of the same quantity in the Yang-Mills version of the theory. 
In fact, for the interaction energy of static quarks in the fundamental repre- 
sentation, we have 

E{xi,X2) = l^i - X2\ 

AnaR 

xi and X2 are the locations of the quarks in one dimension. 
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VI. Partition Function Near Critical Point 

The partition function z in (6) can be written by expanding the cosine 
and the log in their respective power series. The action is 

^ co^O- 1 
-2S = Ej:'-^ + Y^4-^log2[l-cos{e,-9,)] 



N 



^ AT 



\ 0^ df 

1 ^ + -L 

^ AT L 2! 4! 

E]^iog(^^-^.r-E 



+ 

; iV2 3.4 



Kj 

Now performing a scale transformation 

we get 

-2^ = -T7ECV2 + ^Elog|6-0 



N 
+E- 



1 i<j 
A^(A^ - 1) 



2N^ 



6*^"^ stands for the coefficient of e.g., 



'■<3 



12 



^ - N^^' mf^^ 3.4.5.6 
After scaling the measure factor we have 

, = e''"(-4--)i/;jnf n9(V£.-i& 

xe • e ^"=1-=— 
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^=-^Ef + ]^Elog|e.-e.r (33) 

1=1 i<j 

The 9 function in (33) is a cut-off tfiat reflects the compactness of the gauge 
group. 

So far, (33) is exact without approximation. In the large hmit, one 
can once more consider a WKB expansion in 1/iV^. Since this is the same 
problem as before except for a change of scale, we have, for the distribution 
of eigenvalues, 

u{0 = -^u{^/^/E) = 2cos^Vl-^sin2e/2^ (34) 

The function u above is given by (19b), since we are, ultimately, interested 
in large values of E. Expanding, we have 



n=l 



This means that the contribution of the terms involving the operators 9^"'\Cd 
in (33), to the leading term in the WKB expansion in l/N"^, are negligible 
near critical point i.e., as — >■ oo (N.B. the large E limit is taken after the 
large N hmit). Further, the cut-off in (33) is automatically repected, and we 
have 

z{E ^ E*) « e^'(^-^i°^^)i- r n — e^'^ (35) 

We now recall that if if is a iV x iV hermetian matrix with iV^ independent 
components, the integration measure over H is defined by'' 

N 

dH^H dHii n dH^dH^^, Hij ^H^ + iHl^ (36) 

i=l i<j 

Once more one can perform a separation of variables to extract those which 
are left invariant by a unitary transformation. Let W be the unitary matrix 
that diagonalizes H, 

H = W+h W, hij = 5ij^j 

then 

N 

i=l i<j 
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and (35) becomes 



(37) 



Therefore, the partition function of the original gauge theory near critical 
point is 



oc e 



-N^V\ogE/2 



V 



oc 



dH 



-I V 



{VEr-^ ^^^^ 

Expression (38) is what one expects to be proportional to the partition func- 
tion of the Yang- Mills version of the theory, which is defined in a spacetime 
volume V by 



Zym 



_ Tr f dxdtF, 

Jv ' 



2 



(39) 



x.t 



F^y is the field strength; g is the dimensional coupling constant. We have 
fixed the generalized axial gauge 

Ao{x,t) = 

A(a;,to) = 

In this gauge one can express Ai uniquely in terms of the field strength 

Ai{x,t) = Foi{x,t')dt' (40) 

•Ito 

Making this change of variable in (39) 

Zym oa [[dF{x,t)e 

x,t 

At this stage we can manipulate using a cut-off to write 



Zym oc ( / dFe ^a'TiF' 



V 
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dF is the integration measure of the Hermetian matrix F. Introducing the 
dimensionlcss couphng = g'^a'^, which vanishes as a — 0, and the dimen- 
sionless matrix H = Ic^F j \fFi\ E = 2/gQN a fixed number as ^ and 
iV — )■ oo, we have 

Zv« - f (41) 



E) 



Hence the formal equivalence of (38) and (39) is established. 

Using entirely similar arguments it is easy to show that in the thermo- 
dynamic limit of large N, near the critical point, the Wilson loop is given 

by 



W[C] = 



2E\Nt,, 



(42) 



N 



1=1 I 



i i<j i=l 



(43) 



(N.B. to leading order in WKB { / = in agreement with 28). (42) 



can be rewritten as an integral over a random Hermetian matrix H: 

W[C] 



1/1 \ 
1 - — ( -TrH'' 
2E \N 



(44) 



j^TrH^) = I dH(^TiH^) e-f 



/ 



dHe-fT^H^ 



(44) is exactly the expression for the Wilson loop in the Yang-Mills theory. 
The steps are similar to those which led from (39) to (41). 
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VII. Exact Computation of Partition Function 
for Fixed P and 

So far, we have studied the U{N) gauge theory in the hmit iV — )■ oo, 
holding the ratio E — 2/3/N fixed. We now wish to compare these results 
with computations for finite N. 

The partition function z can be expressed as a Toeplitz determinant. 
The trick lies in expressing the repulsive measure factor in (4) as a product 
of Vandermondc determinants 

[](e*^* - e''^^) l[{e-'^' - e"^^^) = detF detF* 

i<j i<j 

Vem = e'^'^ (45) 
dety= ^Ji-JNyiJi-v^j^ 

Jl—Jn 

Substituting (45) in (4) we have 

1 r2n N ^Q, i - ^(-Wi + 2/3 COS Of] 

= detM (46) 
M,,- = /,_,(2/3) = fj ^e^(fc-j)e+2/3cose 

Ini'^P) is a Bessel function with imaginary argument^-*^. 

It is clear that as long as N is finite, z is an analytic function of /3. It 
is also known that the U{N) gauge theory with finite has only 2 trivial 
fixed points of the renormalization group^, (3 = and (3 = oo. As one goes 
to short distances, the coupling constant renormalizes to larger values and 
^ = oo is a critical point near which one can define a continuum theory. 

Prom formula (26) we have, for the Wilson loop, 

W[C] = uj^ 

It is instructive to compute a;(/3, N) near the fixed points /3 = and /3 = oo. 
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A Taylor series expansion of the function 

F^(/3) = N) = det(/,_,(2/3));vxiv (48) 

was done on a computer around the fixed points f3 = and /3 = oo for 
N = 2,3, 4, 5. This calculation was done by Mark Sweeny The result are 

(i) Taylor series around ^ = 0: 

f'3(/3) = l + /J' + y + y + ^/3' + --- (49) 

We see that 

Fn{P) ~ e^' (50) 

is a fairly accurate fit that grows better with increase in N. Prom (50) we 
can compute 

a;(AA^)c^f = f (51) 

(ii) Taylor series around (3 = oo: 
Defining 

The Taylor expansion of /Ar(x) around x = is 

^ / X . n 2 61 3 151 4 771 5 

Ux) = 4x - 7x^ + —x^ - —x^ + —x^ ■ ■ ■ 

fsix) = 128x^ - 912a;^ + 3297x^ (52) 

U{x) = 49152a;^ - 860160x^ + 7452672a;* + ■ ■ • 

Note that the powers of the leading terms are given by N{N — l)/2. This 
leads to 
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2E' 

The small corrections are computable from (52). We note that this result is 
also inferred from the exponential factor that scales z in (33). We do not 
know a simple explanation for the linear fit in (51). 

We see that the finite N results (51) and (53) near the fixed points exactly 
match the corresponding expressions (28) computed in the thermodynamic 
limit of A'^ ^ oo. 

VIII. Conclusion and Discussion 

The analytic computation done in the large N limit holding E = 2f3/N 
fixed is in very good agreement with the computation for finite N, espe- 
cially near the fixed points and it is important to emphasize this agreement 
in the critical region which is the physically relevant one. The difference 
is that while the computation for finite involves analytic interpolation 
between weak and strong coupling the thermodynamic computation shows 
non-analytic behavior ai E — 1. Certainly this non-analyticity is a pathology 
of the N — oo hmit. We do not know whether the quahtative agreement of 
the two methods and even the anomaly ai E — l,is a. feature of the theory in 
higher dimensions which is vastly more complicated. It is desirable to express 
the summation of the planar diagrams in higher dimensions as a coupled set 
of integral equations in the mean field approximation. 

We end with a final comment that the expression of the partition function 
as a determinant implies that this theory has a fermionic representation. 
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